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We investigate the dynamics of a boundary field coupled to a bulk field with a linear coupling in 
an anti-de Sitter bulk spacetime bounded by a Minkowski (Randall-Sundrum) brane. An instability 
criterion for the coupled boundary and bulk system is found. There exists a tachyonic bound state 
when the coupling is above a critical value, determined by the masses of the brane and bulk fields 
and AdS curvature scale. This bound state is normalizable and localised near the brane, and leads to 
a tachonic instability of the system on large scales. Below the critical coupling, there is no tachyonic 
state and no bound state. Instead, we find quasi-normal modes which describe stable oscillations, 
but with a finite decay time. Only if the coupling is tuned to the critical value does there exist 
a massless stable bound state, as in the case of zero coupling for massless fields. We discuss the 
^rr^ t relation to gravitational perturbations in the Randall-Sundrum brane-world. 
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I. INTRODUCTION 



in 

There has been huge interest over recent years in higher-dimensional models involving degrees of freedom which 
propagate on lower dimensional surfaces (or branes) in a higher-dimensional spacetime (or bulk). In particular for 
codimension-one branes in an anti-de Sitter (AdS) spacetime, bulk fields have a zero-mode localised near the brane, 
and hence can yield a lower-dimensional effective theory at low energies without a compact extra dimension [1]. In 

■ traditional Kaluza-Klein compactification the low energy effective theory arises from the lowest eigenmodes of the 
, discrete spectrum of bulk fields. By contrast the lower-dimensional field living on the boundary of a bulk spacetime 

■ can in principle couple to an infinite tower of bulk states even at linear order in perturbation theory [2] . 

In practice the more massive bulk modes may be weakly coupled to the boundary at low energies and hence the 
*n . tower may be neglected. But if we want to determine the quantum vacuum state of the boundary field at all energies 
we need to be able to deal with coupling to bulk modes. This is a practical problem in cosmological models of 
<i , inflation driven by an inflaton field on a brane [3] . The large-scale structure of our universe is supposed to originate 
"^T* ' from small-scale vacuum fluctuations of the inflaton field during inflation. But these inflaton fluctuations on the 
Oh, brane start on arbitrarily small scales where they are strongly coupled to the infinite tower of bulk graviton modes 

■ [4]. Therefore, a truncated tower may not give a useful description of the vacuum state on small-scales from which 
the large-scale structure of our universe is supposed to be derived. 

To derive the vacuum state of a coupled boundary-bulk system we must include the full spectrum of normalisable 
modes of the system. We note that the spectrum of normalisable modes can include discrete bulk eigenmodes with 
complex eigenvalues in addition to the familiar modes with real eigenvalues that form a complete orthonormal set 
of basis functions for any classical solution. Dubovsky et al. [5] studied quasi-local states with complex eigenvalues 
for massive bulk fields in the anti-de Sitter bulk bounded by a vacuum Randall-Sundrum brane-world. This work 
has been extended to include the curvature of the brane [6] , and the self-coupling of the bulk field on the brane [7] . 
Recently, Seahra [8] has identified the quasi-normal modes of gravitational waves by seeking complex eigenvalues of 
the bulk mode equation. As far as we are aware there has be no attempt previously to identify the bound states or 
quasi-normal modes of coupled boundary-bulk field theories in an anti-de Sitter spacetime. (See Ref. [9] for a different 
approach.) 

Cartier and Durrer [10] did consider the general solution for gravitational wave modes coupled to an anisotropic 
stress on the brane in a Randall-Sundrum model. They found that the system could describe normalisable modes 
with the anisotropic stress growing exponentially with respect to time on the brane. They went on to suggest that 
instabilities might be generic in any orbifold model for the bulk. But the generality of their argument is open to doubt 
without including the self-consistent evolution of a reasonable matter source on the brane, and without checking that 
the instability is not driven by energy flux from the past Cauchy horizon. 

As a first attempt to study the dynamics of a boundary field coupled to bulk fields in AdS we study a simple model 
of massive bulk and boundary fields with a linear interaction term on the brane defined in Section II. We begin 
by reviewing in Section III the recent work of George [11] who considered an oscillator on the Minkowski boundary 
coupled to a massive field in a Minkowski bulk. We then go on to extend this to a massive field on a Minkowksi brane 
in an AdS bulk in Section IV, identifying resonant modes corresponding to bound states and quasi-normal modes, 
and present a condition for the existence of tachyonic instabilities. We explore analytic and numerical solutions for 
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these resonant modes in limiting cases in Section V. We summarise our results in Section VI and discuss possible 
implications for gravitational perturbations in Randall-Sundrum brane- worlds. 

II. THE MODEL 

In the (d + l)-dimcnsional bulk spacctimc V, with metric Gmn, we consider a free scalar field <f>, with Lagrangian 

= ~\g mn 4>, M <t>,N - \m 2 ^ 2 , (2.1) 
and on the d-dimensional boundary <9V, with metric g^, we consider a field q with Lagrangian 

C q = - l -g^q 4l q. v - l -^q 2 . (2.2) 
For simplicity we consider a linear interaction between the boundary field and the bulk field on the boundary 

Ant = ~PH ■ (2.3) 

Note that the system is symmetric under (3 — > —(3 and <j> — > —<f> or q — ► —q. Henceforth we take (3 > without loss of 
generality. The total action is thus 

S= f d d+1 xV^G£t + f d d x^g~[C q + £ int ] . (2.4) 

JV JdV 

The coupled wave equations are then 

Wnq = ij, 2 q + f3<j> , (2.5) 
on the brane, where <po denotes the value of <f> at the boundary, and a free wave equation for <j) in the bulk 

( d+ Vu4> = m 2 <}>, (2.6) 

subject to the boundary condition [17] 

$ = f?, (2-7) 

at the brane, where <p' signifies the normal derivative at the boundary. 

Note that if we were considering a linear interaction between two massive boundary fields with total Lagrangian 

C = £ qi +£ q2 -Pq iq2 , (2.8) 
then we could easily diagonalisc the mass-matrix and identify the effective squarcd-masscs 



2 ml + ml ± \/{m\ + nii,) 2 + 4/3 2 - km\m\ 

m± = - , (2-9) 

on of which becomes negative for 

f3 2 > m\m\ , (2.10) 
signalling a tachyonic instability. Thus, we should not be surprised to find an instability for large coupling. 

III. MINKOWSKI BULK 

George [11] recently considered an oscillator on a flat boundary coupled to a field in a flat two-dimensional Minkowski 
spacetime. We write the corresponding bulk metric in (d + l)-dimensional Minkowski spacetime as 

ds 2 = -dt 2 + d'',l.r,d.r , + dy 2 . (3.1) 

The coupled wave equations (2.5) and (2.6) thus become 

q + k 2 q = - A i 2 <Z-/#0, (3.2) 

4> + k 2 (f> = (f>"-m 2 (/), (3.3) 

where a prime denotes a derivative with respect to y and a dot with respect to t, subject to the boundary condition 
(2.7), and we consider a single Fourier mode with wavenumbcr k on the (d — l)-dimcnsional sub-space. 
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A. General solution 

Formally the general solution can be written as a sum over modes [11] 

/oo 
d P C p T p (t), 
-OO 

/OO 
dpA p e^T p (t). 
-OO 

The bulk wave equation (3.3) requires 



T n = -uiT, 



p > 



where, uj p is given by the dispersion relation 



lo 2 = p 2 + m 2 + k 2 



(3.4) 



(3.5) 



(3.6) 



In addition the boundary wave equation (3.2) and the boundary condition for the bulk field (2.7) require that the 
following consistency relations between the coefficients 



(p 2 +m 2 - p?) (C p + C_ p ) = p(A p + A_ p ) , 
ip(A p -A. p ) = l{C p + C^ p ). 



(3.7) 
(3.8) 



Eliminating (C p + C_ p ) from these two constraints we obtain a relation between the coefficients in the bulk mode 



■ l 2 i 2 2\ 
ip{p + rn — p, ) 



2 



-4, 



ip(p 2 + m 2 - p 2 ) + 



d 2 



A_, 



(3.9) 



The constraints interpolate between Neumann boundary conditions at weak coupling (/3 — > 0) and Dirichlet boundary 
conditions at strong coupling (J5 2 — > oo) [11]. 



B. Resonant modes 



Conventionally, p is assumed to be real; however, resonant modes with complex p play a crucial role in a scattering 
process. Resonant modes can be found by demanding that the bulk field is purely out-going at infinity y — > oo. We 
take the temporal mode function as 

T p (t) = e-^"\ (3.10) 
and impose the condition cj)(y,t) — > e ?py e _iWpt as y — > oo. In order to satisfy this condition, we must have 



A- p = 0. 



(3.H) 



Substituting this constraint into equation (3.9) yields 



p{p 2 +m 2 - p z )+i'— =0. 
The out-going boundary condition are satisfied when the real parts of p and uj p have opposite signs, i.e., when 



(3.12) 



5R(p)5R(w p ) > 0. 



(3.13) 



Resonant modes correspond to the poles of the Green function [5]. Let us compute the 5D Green function 
A(x fJ -,y 7 x fJ - ,y') that satisfies 



(9y - p 2 ) A (p, y, y') = S(y - y') , 



where we performed a Fourier transformation along the brane 

A(^,y;^',y') = 



;d ^_ e i P .(x-x') A(/9)2/)y0 . 



(27T) 



(3.14) 



(3.15) 
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and p 2 



-P 



uj„ — m 



k 2 . The boundary condition for the Green function at the brane y = should be 



determined by the boundary conditions for the scalar field in question: 

4>'{p) = f 

<f>(p) = ^(p 2 + m 2 - p 2 )q(p), 



(3.16) 



where <f)(p) and q(p) are the Fourier modes of </>(0, x^) and q{x^) respectively. In order to satisfy these boundary 
conditions, we must impose the boundary condition for the Green function as 



(p 2 +m 2 -p 2 )d y A-^-A 



= 0. 



(3.17) 



y=o 



We also impose the out-going boundary condition at infinity. The Green function has a simple form if one of the 
arguments is on the brane; 



A(p,y,0) 



2 i 2 2 
P + TO — p 



aiPV 



ip(p 2 + m 2 — p 2 ) — *y 



(3.18) 



We notice that the condition for Eq. (3.12) is the condition for the poles of the Green function. 

The causal boundary condition for the Green function is determined by the way in which the contour of integration 
is closed around the poles. In order to impose the retarded boundary condition we must close the integration contour 
in the upper half complex uj p plane for t < t' and subtract the contributions of the poles in the upper half complex 
ujp plane. For t > t', we can close the integration contour on the lower half of the complex lu p plane and the late 
time behavior of the scalar field can be understood by investigating the structure of singularities such as poles in the 
integrand. 

Fig. 1 shows the solutions for lo p in a complex uj p plane for p > to and /3 > 2mp 2 . The integration contour must 
be taken like C in order to impose retarded boundary condition. We find two kinds of poles described below. 



Im( (Op ) 



c 



Re( (Op ) 



FIG. 1: Location of poles in a complex ui p plane for \i > m and (3 > 2mp 2 with k = 0. Filled circle show poles with the 
out-going boundary conditon and open circles show poles with the in-going boundary condition. The contour is chosen so that 
the retarded boundary condition is satisfied in the case of the out-going boundary condition. 



1. Bound state 



The bound state is a mode that is exponentially decaying and localised near the brane. Setting p = ia we see that 
we are seeking roots of the cubic equation 

B 2 

a 3 + cr(p 2 - to 2 ) - '— = 0. (3.19) 
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with positive real part. 

Equation (3.19) always has one (and only one) positive real root [11]. This will describe stable oscillations if u> 2 > 
which, from Eq. (3.6), requires a 2 < m 2 as k — -> 0, or equivalcntly [11] 

[3 2 < 2m/i 2 . (3.20) 

On the other hand the bound state describes an unstable mode which can grow exponentially in time above a critical 
coupling /3 2 rit = 2m/i 2 . In this case, we have two poles on the imaginary axis in the complex iv p plane: a pole in the 
upper half plane describes an unstable mode. 



2. Quasi-normal modes 



In addition to the bound state, there are solutions with complex uj p located in the lower half lu p plane. Thus, 
these solutions lead to an exponentially decaying mode in time. We should note that this conclusion depends on the 
boundary condition. If we impose an in-going boundary condition so that <f>(y,t) — > e~ lpy , e ~ llAj p t with the conditions 
Eq. (3.13), poles are located at the upper-half complex lu p plane, leading to an exponentially growing mode in time. 
However, a physical boundary condition, i.e., the out-going boundary condition under which there is no energy flux 
from infinity, can eliminate these unstable modes, therefore, these instabilities do not describe a physical instability 
of the system. On the other hand, an unstable bound state for j3 > /3 C rit cannot be eliminated by a choice of the 
boundary condition, indicating a physical instability of the system, i.e., a tachyonic instability. 



IV. ANTI-DE SITTER BULK 



We now extend the analysis to consider a boundary field coupled to a massive bulk field in an Anti-dc Sitter (AdS) 
bulk spacetime, as in the Randall- Sundrum brane- world [1]. We write the corresponding bulk metric, with Anti-de 
Sitter curvature scale £, as 

I 2 

ds 2 d+1 = e - 2M/e (-dt 2 + fidxidxj) +dy 2 = -j (-dt 2 + S i: >dx l dx J + dz 2 ) . (4.1) 

Note that the boundary is now located at z = i, and the surface I < z < +oo at fixed time t defines a Cauchy surface, 
with t — > — oo the past Cauchy horizon [18]. 

The coupled wave equations (2.5) and (2.6) become 

q + k 2 q = -,! 2 g -/# , (4.2) 
3 m 2 f 2 

t + k 2 = <t>" -U' -^i-<t>> (4.3) 
z z z 

subject to the same boundary condition (2.7) as before, where a prime now denotes derivatives with respect to z, and 
we again consider a single Fourier mode with wavenumber k on the (d — l)-dimensional sub-space. 



A. General solution 



Formally the general solution can be written as a sum over modes 



«(*) 



dpCpT p (t) 



= / dp [A p HP(pz) + B p Hl 2 \pz)} T p (t) 



(4.4) 



where and are Hankel functions of first and second kind and order v = \/A + m 2 £ 2 . The bulk wave equation 
(4.3) then yields the time dependence for T p (t) 



-u}T~ p , 



(4.5) 
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where is given by the dispersion relation 



uj = (? + k 2 



(4.6) 



The boundary wave equation (4.2) and the boundary condition for the bulk field at the branc (2.7) impose the 
conditions 



(p 2 - p 2 )C p = /3 \ApHf\pl) + BpH (2 \pl) 



£c~ - A- 

2 O p — Zip 



^H£\p€) + pHl\{p€) 



^Hi 2 \pi)+pH^\(pi) 



(2) 



(4.7) 
(4.8) 



from which wc eliminate Cp to get 



-Bp 



Kf-^H^W) + { 2 -j^(p 2 -p 2 ) + ^}H^(pi) 
p(p 2 -p 2 )Hl 2 \(p£) + {^(f?-p?) + t} H ?\pl) 



(4.9) 



We have used various well-known relations for Bessel functions and their derivatives which can be found in Ref. [12]. 



B. Resonant modes 



In order to find resonant modes, we impose that the modes become plane waves at large z, 
where 9 is a phase and Kp is a constant. Then we must have Bp = and Eq.(4.9) gives 



p(p'- t i 2 )Hi 1 \(pi) + 



2-v 



2 1 



HW(pe) = o, 



(4.10) 



(4.11) 



with 5R(/5) > 0. Note that if we had instead imposed the condition would have obtained the 

condition 



p( P 2 -p 2 )Hl 2 } 1 (pe) + 



2-v 



I 2 



HW(p£) = 0. 



(4.12) 



The solutions to which are just the complex conjugates of the solutions to Eq. (4.11). 

Again, solutions to Eq. (4.11) are related to poles of the Green function. Let us compute the 5D Green function 
A(x M ,z;a; A ' , z') that satisfies 



where we performed a Fourier transformation along the brane 

d 4 p 



f! (dl _ l da /j A (p, z, z>) = (f ) * *(* " z>). 

e ip-(x-x') A (fr ZjZ >^ 



A(^,z;x^ ,z') 



(2tt) 4 



(4.13) 



(4.14) 



and p 2 = —p 2 = uPp — k 2 (compare to Ref. [13]). The boundary condition for the Green function at the brane z = I 
is given by 



(p 2 -p 2 )d z A-^-A 



0. 



(4.15) 



We also impose the boundary condition at the Cauchy horizon of the AdS spacetime so that the waves are out-going. 
The Green function has a simple form if one of the arguments is on the branc; 



A(p,z,£)=(£) ; 



(p 2 - p 2 )H^(pz) 



(4.16) 



p(? - M 2) ^ + ^ ^ + £j H (D ^ 

As in the Minkowski case, the condition for the bound states is the condition for the poles of the Green function 

AQW)- 1 ^. (4.17) 
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1. Bound states 

The AdS curvature in the bulk dramatically changes the condition for the existence of a bound state. Setting 
p = id, we seek the solution for Eq. (4.11). We can use the relation 

K v {x) = l ^e i ™l 2 HV(ix) (4.18) 

where K v {x) is a modified Bcsscl function [12], to write Eq. (4.17) as 

Ku-^ai) (u-2) (3 2 1 



K v (al) I 2 (a 2 +p 2 ) 



0. (4.19) 



The existence of a positive real root depends on the value of the coupling. There is always one (and only one) positive 
real root for f3 > /3 C ritj where 

Pl it = 2{v - 2) ( i 2 r 1 . (4.20) 

As in the Minkowski bulk case, this describes an unstable mode with w| < when k = 0. Note that /3 cr it becomes 
/3 cr it = 2m 2 p for ml 3> 1, which agrees with the critical coupling in the Minkowski case. For j3 = /3 CT it, there is 
a massless bound state a = 0. Unlike the Minkowski bulk case, there is no positive real root for a small coupling 
P < Aadti instead we find quasi- normal modes with complex p. 

2. Quasi-normal modes 

We have a rich structure of quasi-normal modes for an AdS bulk. Indeed, the bulk scalar field has an infinite number 
of quasi- normal modes even without a coupling to the brane field [8] . As was shown in Ref. [5] , if we introduce a mass 
to bulk fields in AdS spacetime, we find quasi-localized states with finite width. Quasi-normal modes with complex 
eigenvalue p describe metastable states that decay into continuum modes within a finite time. With zero coupling, 
[3 = 0, the condition for quasi- normal modes for the bulk scalar field is given by 

H^\(pt) 2-v , 

=°- (4 - 21) 

At low energies ml <§C 1, the solution for p is given by 

2 

^ = HL +iT} r = -^m 2 (m^) 2 , (4.22) 
z lb 

with 3?(p) > 0. If we take the temporal mode function to be Tp(t) = e ~ lw p*, we must impose 3?(cJp) > to satisfy the 
out-going boundary condition. On the other hand, if 5i(wp) < 0, we must impose <f> oc z z / 2 e~ % P z with 3?(/5) > 0. A 
solution for p at low energies becomes 

2 

P 2 = \ - *T, r = -^m 2 (m£) 2 , (4.23) 

with SR(p) < 0. 

These poles are located in the lower half cu p plane (see Fig. 2) and these solutions lead to exponentially decaying 
modes in time. Note that we could have exponentially growing modes in time if we imposed an in-going boundary 
condition, just as in the Minkowski bulk case. However, an in-going boundary condition implies energy flux from the 
past Cauchy horizon of the AdS spacetime. Therefore, this is not a physical boundary condition. 

A non-zero coupling to the brane field changes the location of the poles in a complicated way. We will study the 
effect of the coupling in the next section in various cases. 
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FIG. 2: Location of poles in a complex u p plane for massive scalar fields with k = 0. Filled circles show poles with the out-going 
boundary condiiton and open circles show poles with the in-going boundary condition. There are an infinite number of poles. 



V. BOUND STATES AND QUASI-NORMAL MODES IN LIMITING CASES 



A. High-energy limit 



Let us first consider the case where ml 3> 3> 1 and search for solutions with pi 1. In this limit, we can 
show 

—^— + — ^i^p—^. (5.1) 

The condition for the bound state becomes the same as Minkowski bulk case, equation (3.12) if we make the identifi- 
cation 

p 2 = p 2 ~m 2 . (5.2) 

Thus, we recover the Minkowski results in the limit when the AdS length scale is much larger than the Compton 
wavelengths of both the boundary and bulk fields. 



B. Low energy limit 

Let us now consider the low-energy limit where mass scales are much less that the AdS scale (ml <C 1 and p£ <C 1) 
and search for bound states with \pl\ <§; 1. We can expand the Hankel functions for small arguments and the condition 
for a bound state (4.11) becomes 

(P 2 P 2 ) (V + in^f ~ ! y)= ^0*, (5-3) 

where we assumed m ^ 0. At weak coupling. 1 ji — > 0, there are two independent fields. One is the brane field with 
p = p and the other is the quasi-bound state of the bulk field with Eqs. (4.22), (4.23) 

Neglecting the finite decay width, suppressed by ml <C 1, we can construct the effective coupled equations for the 
brane field q(t) and the bulk fields on the brane <p(Q, t) 

4> + k 2 cf> + ^m 2 (f> = -r 1 ^, 

q + k 2 q + p 2 q = -j3(j). (5.4) 

This d-dimensional effective theory coincides with the example of two massive boundary fields with linear coupling 
defined in Eq. (2.8) if we identify m\ = p 2 , m\ = m 2 /2 and /? = (3/1. We would expect to be unstable above a 
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critical coupling given by Eq. (2.10) which corresponds to 1 > £m 2 p 2 /2. Note that this critical coupling agrees with 
Eq. (4.20) for mi < 1. 

However, this effective theory neglects the slow decay indicated by the imaginary part of the effective mass p 
suppressed by pi <C 1. In the following sections, we study the effect of this small imaginary part induced from the 
curvature of the AdS spacetime. 



C. Massless bulk and brane fields 

Let us consider the simple case where to = /i = 0. Henceforth, we will use units where i = 1, so that the consistency 
relation, Eq. (4.11), for massless fields can be written as 

^i (1) (p)-yf(p) = 0. (5.5) 

When /3 = there is no coupling and there is a zero mode solution p = 0. If the coupling is turned on, this will 
no longer necessarily be a solution. For small f3 there should be a root near p = 0, so we will investigate how the 
consistency relation is perturbed. The Hankcl functions can be expanded for \p\ <C 1 by 

rrfil/M 4i i i ~ f p\ 2ir + 4Ci — 3i „ 2 „. 

H ^ = -w~-^-* p l0 H2j + — is? — p " + -"' (5 - 6) 

r( 1 )/"\ 2 ^ i * -i ( P\ n + 2Ci — i 



H\ L, {P) = —,+-p\o g r-)+— p + ..., (5.7) 

irp tt \2 J 2n 

where C « 0.577 is the Euler constant [12]. So the consistency relation can be approximated as 

2 *~2 i -4, /'A 7r + 2Ci-i„ 4 2i/? 2 i^ 2 i o2 „ 2 , /A 3i - 4C*i - 2tt „ 2 . . 



tt r 7r r ""°V 2 y 27T r ' Trp 2 ' 2tt 87r'" r ""°V 2 y ' 32tt 

To lowest order in p this has solution 



P i = [3 2 => P 2 =±P => / 5 = e" n/4 v //3, n = 0,1,2,3. (5.9) 
if we continue the solution p = \/]3 by adding higher order terms in (3 we get 

P = + ^/3 3/2 log/? + i(C - log2)/3 3 / 2 - ^/3 3 / 2 + • • • . (5.10) 
Using the dispersion relation (4.6) we see that the frequency is not real, but is given approximately by 



/3 



2 



u = ±JeTp---g== (5.11) 

again, keeping only the lowest order real and imaginary parts. Note that for Sft(wp) < 0, we must use the solution for 
Eq. (4.12) to impose the out-going boundary condition. We see that the poles are located at a lower-half plane and 
these solutions are quasi-normal modes that describe massive metastable states. For the root near iy/]3, we get 

P = iVTl-*-^P 3/2 log/3 - i\(C - lo g 2)/? 3 / 2 , (5.12) 

where the real part vanishes. This is the massive bound state which leads to a tachyonic instability. 

The consistency relation (4.11) can be solved numerically for any given values of the parameters to, p and /?, 
normalized by setting I = 1. The roots for massless fields and small coupling, (3 = 0.01, can be seen on in Fig. 3. 
The right hand plot shows the roots near the origin already calculated analytically. There are also other roots, one 
at —0.419 — 0.577z and a sequence with almost identical imaginary parts of — 0.355i exactly analogous to those found 
by Seahra [8]. 
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FIG. 3: These plots, generated by the Maple 9.5 software package, show as poles the roots of the consistency relation when 
p = m = and (3 — 0.01. They are contour plots of the logarithm of the absolute value of the l.h.s. of Eq. (4.11). The roots 
near the origin can be determined by Maple's root finding routine as ±(0.0997 — 1.40 x 10 _4 i) and ±0.100i. 



D. Massive brane field 



Let us now consider the case were p ^> 1 where we are still considering j3 <C 1. When (3 = there arc roots to 
Eq. (4.11) at p — ±/x. When [3 ^ we can work out how these are perturbed. By writing p = p + 5p we can expand 
the consistency relation by a Taylor series about p = p to get the lowest order correction 

S P = 7 77\ 77\ ~ ( 5 - 13 ) 

Ap[pH i l }l l (p) + (2- V )H^{p) 

So the roots are 

± ,; + r&w | , (5 . 14) 

<pU™,Wt (2 -v)nP(p)) 



and the frequency u> will be given by the dispersion relation (4.6). 
One case of special interest is when m = 0, where this simplifies to 

^ ± L + £,^M). (5.15) 

The imaginary part always has the opposite sign to the real part. We can use the well known asymptotic behaviour 
of Bcsscl functions [12] to write 



'h£\p)\ r-f/i p^o 



so for a light field on the brane, where j3 <C p <C 1 



(5.16) 



>"*("-£')■ "- ± ^ f -i7inT? < - (5 - 17) 

There is also a solution with purely imaginary p which corresponds to the tachyonic bound state. 



11 



VI. CONCLUSIONS 

In this paper we have investigated the behaviour of a boundary field theory on a Minkowski brane linearly coupled 
to a bulk field in anti-de Sitter spacctimc. We have calculated an instability criterion for the coupled boundary and 
bulk oscillators. Above a critical coupling CT - lt = \/2(v — 2)p 2 /£ there is a state with a purely imaginary effective 
mass. This state is spatially bounded and normalizable in the bulk, and describes a tachyonic instability on large- 
scalcs (small k). When either field is massless /3 cr it = 0, so that the instability will always occur for non-zero coupling. 
Below the critical coupling there is no tachyonic state and no bound state, in contrast to the case of a Minkowski 
bulk [11]. Instead we find quasi-normal modes. 

The quasi-normal modes describe the late-time behaviour of the system when there is no tachyonic instability. For 
weak coupling, there is a quasi-normal mode which represents a very slowly decaying oscillatory solution. For instance 
in the low energy limit (ml <C 1) we find a decay time of the order (mt)~ 1 m^ 1 . This contradicts the conjecture of 
Cartier and Durrer [10] that orbifold boundary theories necessarily have instabilities. It is necessary to impose a 
physical boundary condition in order to determine which poles in the complex frequency plane describe the behaviour 
of the system. Quasi-normal modes are determined by a purely outgoing condition at infinity and thus only poles in 
the lower half complex plane describe the late-time evolution of the system. 

Quasi-normal modes are not normalizable modes. For the classical problem normalizability is not a concern because 
the quasi-normal mode solution is only valid in the causal future development of the brane (see Ref. [14] for a 
dicussion). On the other hand normalisability is a concern if one wishes to construct an initial quantum vacuum in 
terms of independent oscillators with finite action. In contrast to the work of George in a Minkowski bulk [11] we are 
not able to give the full spectrum of the quantum vacuum state below the critical coupling. 

When the coupling is zero, a massless bulk field behaves in a way very similar to gravitational waves in the 
Randall-Sundrum brane- world [1]. There is a massless zero mode solution on the brane and a spectrum of massive 
Kaluza-Klein states. In addition, there is a spectrum of quasinormal modes found by Seahra [8] which describe the 
evolution of the gravity waves at intermediate times on the brane, though the zero mode dominates at late times. 
When a non-zero coupling is introduced this behaviour is modified dramatically. The zero mode, p = 0, is no longer 
a solution but is perturbed to solutions with complex p. Without the coupling, the boundary condition at the brane 
is purely Neumann, thereby admitting a bounded zero-mode. One can think of the coupling as mixing Dirichlet and 
Neumann boundary conditions. A similar effect could occur with a moving boundary, such as an expanding brane 
universe. 

An important difference in the case of gravitational perturbations in the bulk is that these are not in general coupled 
to linear matter perturbations on the brane and metric backreaction is treated as a second-order effect. An important 
exception is the case of fluctuations in an inflaton field driving inflation on a brane. Inflaton fluctuations on the brane 
are linearly coupled to bulk metric perturbations. These fluctuations start on small scales where they are strongly 
coupled to bulk gravitational perturbations [4]. Thus, we must specify the vacuum state for a coupled boundary 
(inflaton fluctuations) and bulk (metric perturbations) system [15]. 

In order to address this issue, we need to extend our analysis to the case of a de Sitter brane. Due to the curvature 
of a de Sitter brane, there appears a mass gap [16] in the spectrum of the free bulk field and there arises the possibility 
of having stable massive bound states in the mass gap [7] . One could crudely model the existence of a cosmological 
horizon in de Sitter space by imposing an infra-red cut-off on wavelengths larger than the Hubble length, suppressing 
the tachyonic instability for sufficiently small coupling or large Hubble rate. We leave a full analysis of the de Sitter 
brane for future work but it is intriguing to note that a de Sitter brane might provide a good model for a coupled 
branc-bulk quantum field theory, where a Minkowski brane does not. 
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